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ABSTRACT 



Synthetic Aperture Radar wind retrieval models have great potential to 
accurately depict the mesoscale wind field on the order of hundreds of meters. 
However, there are still significant hurdles to overcome in applying the theory in 
a fully automated algorithm. This thesis is a result of collaboration between the 
Nansen Environmental and Remote Sensing Center, Bergen, Norway and the 
Naval Postgraduate School. It examines two wind retrieval models on basis of 
in-situ and remote data from NORCSEX ‘95 in the interest of improving our 
understanding and application of these models. The individual models are based 
on Bragg backscatter from the ocean surface however, the CMOD4 model is 
directly related to the backscatter while the SWA model is related to the spectral 
resolution of the backscatter field. Each model has specific advantages and 
disadvantages related to processing and automation. For this data set they show 
very good agreement with each other and with in-situ measurements. Plots of 
optimum wind vectors derived from a combination of CMOD4 and SWA wind 
vectors are shown to illustrate the high resolution wind fields available with this 
technology. 
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L INTRODUCTION 



The capability to sense remotely the surface wind field over the ocean at a resolution 
on the order of hundreds of meters would be of great value to Naval and civilian 
meteorologists, especially in coastal areas. Scatterometers employed on spacecraft are limited 
in resolution to about 50 km. It is not practical to improve this resolution since it would 
require a much larger antenna. This physical limitation on antennae size can be relaxed by 
creating a synthetic aperture through analysis of the frequency shift of the transmitted signal. 
SEAS AT in 1978 carried the first operational spacecraft-based Synthetic Aperture Radar 
(SAR), which was able to resolve images of about 25 m . The goal of the SEASAT SAR was 
to obtain high resolution images of the sea surface and sea ice. Operationally, it had limited 
coverage in the Northern Hemisphere and no coverage in the Southern Hemisphere because 
there was no on-board storage for the high data rate SAR which required direct real time 
transmission to the four ground stations equipped to receive it. The SEASAT SAR far 
exceeded its goals for providing high resolution images of ocean waves, internal waves, 
underwater topography, eddies, fronts, natural slicks and oil slicks (Beal et al. 1981, Allan 
ed.1983). 

Analysis of the surface wind was not one of the original/primary goals of the SEASAT 
SAR but one that evolved once the images became available. Wind fronts were clearly evident 
in many coastal SAR images where the scatterometer showed a smooth field or was degraded 
due to interference from the land backscatter (Allan ed. 1983). Subsequent SAR’s have been 
employed on the Space Shuttle, ERS-1, ERS-2, and RADARSAT and have been used to 
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analyze the surface wind. However, the exact relationship between the SAR return and the 
ocean wave field is still not completely understood leading to errors in the processed wind 
fields. 

The Norwegian Coastal Shelf Experiments (NORCSEX), located off the southwest 
coast of Norway, have been conducted since 1991 in support of the European Space 
Agency’s European Remote Sensing (ERS) program. The most recent of these, NORCSEX 
‘95, was conducted during the ERS-2 calibration trails. It collected in-situ data from the 
research vessel (R/V) Hakon Mosby, and from an instumented buoy while both the ERS-1 
and ERS-2 SAR’s imaged the same area. This thesis is a result of collabortion between the 
Nansen Environmental and Remote Sensing Center, Bergen, Norway, and the Naval 
Postgraduate school. The Naval Postgraduate School’s participation was to obtain high 
quality in-situ wind data for comparison. It examines two SAR wind models vs. the in-situ 
wind field with the goal of improving our knowledge of the errors in wind retrieval 
techniques. First, background descriptions of SAR and NORCSEX ‘95 are given in Chapter 
II. Descriptions of the wind retrieval models are presented in Chapter III. The results are 
presented in Chapter IV, and conclusions and recommendations for further study in Chapter 
V completes the thesis. 
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II. BACKGROUND 



A. SYNTHETIC APERTURE RADAR 

1. Principles 

The principles behind the operation of synthetic aperture radar are basically the same 
as those of any imaging radar with the exception of the synthetic aperture of course. It is not 
within the scope of this thesis to go into all the details of SAR operation and processing. 
However, the operation of a standard imaging radar will be reviewed and then the formation 
of the synthetic aperture and problems associated with SAR interpretation will be discussed. 

The typical radar uses a pulse of microwave radiation to illuminate an object and 
measures the strength of the returned signal and the elapsed time, which are then processed 
and recorded for display and analysis. The factors affecting the image resolution are the 
pulse width, frequency, range to target, and the aperture or antenna size. When the radar is 
moving with respect to the target, as is the case with satellite imaging radars, then the pulse 
width determines the resolution in the range (cross track) direction and the aperture size 
determines the resolution in the azimuth (along track) direction. A typical geometry for a 
satellite based SAR is shown below in Figure 2. 1 . 

Any radar using the same aperture for both transmitting and receiving has an angular 
resolution (<j>) which is dependent upon it’s aperture length (L) expressed in wavelengths. 

4> = M2L W 



3 



To get higher resolution the aperture length is increased. Present scatterometer resolution 
is limited by its aperture (antenna) size to about 50 km. Doubling of the aperture (antenna) 
size of a satellite borne radar would be very expensive and would result in only a 25 km 
resolution. This physical limitation is overcome by the synthetic aperture radar which uses the 
spacecraft motion along an arc to create a large “synthetic aperture”. The design image 
resolution of 30 m for the ERS-2 is achieved by transmitting a series of coherent pulses while 
the spacecraft travels the distance of the synthetic aperture. These pulses are then corrected 
for Doppler shift and processed to appear as one signal. The time it takes the spacecraft to 
travel the length of the synthetic aperture is only about 0.5 seconds. Fortunately this is much 
shorter than the period of the short gravity waves which SAR images so the error induced by 
this processing method is minimal. 

In the case of a satellite-based SAR, such as ERS-2, the relationship between angular 
resolution (<]>), slant range (h), SAR wavelength (X), azimuthal resolution (r a ), and aperture 
length is: 



r a = $*h!2L (2) 

Much effort has gone into the interpretation of spacecraft SAR images since the first 
SEASAT images in 1978. It is generally accepted that the return signal over the ocean 
surface is due to Bragg scattering. The Bragg model states that the dominant backscattered 
energy will arise from the surface spectral components that resonate with the incident wave. 
The ocean surface can be approximated by a series of planar facets each tangential to the 
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actual surface. The small scale roughness is superimposed on this on this surface so that the 
incident wave scatter has components due to the local slope (facet scattering), and the 
roughness (Bragg scattering). (Curlander and Mcdonough, 1991). The resultant backscatter 
curve is a function of incidence angle (Figure 2.2). From the Bragg theory, the optimal 
wavelength for imaging ocean waves is defined by the relationship: 

A/2 = A sin 0, (3) 

where A r is the wavelength of the Bragg resonant wave. Using ERS-2 data this is 
approximately 5 cm which means that the waves being observed are the small surface gravity 
waves. Experiments have shown the strong correlation between the SAR return and the 
ocean wave field as well as the surface wind field for wind speeds greater than 3 m/s (Barrick 
and Peake 1968, Napolitano 1991 ) 

2. Importance of Ocean Long Waves on SAR Data 
It is believed that the ability of the SAR to detect ocean waves longer than the SAR 
wavelength depends on the modulation of the Bragg scatterers by these waves, although the 
exact relationship is not yet known. Since one of the wind models examined later is related 
to ocean long waves, a short review is useful. There are three principle mechanisms by which 
this modulation occurs: surface tilting, velocity bunching and hydrodynamic interactions 
(Alpers, Ross, andRufenach, 1981). 

a. Surface Tilting 

Surface tilting results from the change in geometric presentation of the Bragg 
scatterers to the satellite. As depicted in Figure 2.3, the swell on which the Bragg waves are 
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superimposed changes the incidence angle and results in the tilt modulation phenomena. The 
radar backscatter is largest on the slope of the long swell wave directed toward the radar. 

The magnitude of the tilt modulation transfer function (MTF) is dependent on both 
the viewing angle of the radar and the incidence angle of the long swell waves. It has been 
demonstrated that the range traveling long waves (waves traveling normal to the satellite) are 
likely to have an MTF about an order of magnitude larger than the azimuthally traveling 
waves (waves traveling in the same direction as the satellite) as seen in Figure 2.4 (Alpers et. 
al 1981). 

b. Hydrodynamic Interactions 

The hydrodynamic interactions are complex interactions between waves and 
currents which result in increased short wave amplitudes and wavelengths near the crests of 
long waves (Figure 2.3). Although it is not well understood, Alpers and Hasselman (1978) 
used conservation of density in the wave interactions to develop a MTF which can be used 
for SAR processing. They also showed that the modulation due to hydrodynamic interactions 
was about half the magnitude of that due to tilt and that its maximum value was also for 
range traveling waves. 

c. Velocity Bunching 

Velocity bunching is the mechanism by which SAR can resolve waves that a 
standard imaging radar of the same resolution cannot. This is because the SAR uses Doppler 
shift to define the azimuth position of a target and any target motion will shift this position. 
The scattering elements on the ocean surface have a time varying radial velocity and their 
non-uniform azimuthal displacement results in wave-like patterns in the SAR images. This 
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mechanism also explains the appearance that ships are azimuthally offset from their tracks. 
If a wave is close to the range direction, the velocity bunching effects are small and linear. 
For waves traveling in the azimuthal direction, velocity bunching is the dominant modulation 
effect. There is only a small angular interval around the range direction where velocity 
bunching is a linear process and similar in magnitude to the tilt modulation (Figure 2.5). 
Alpers et al. (1981) discusses the three modulation mechanisms in detail and develops an 
MTF for the three that is valid in the fairly restricted conditions under which they can all be 
treated linearly. 

The velocity bunching mechanism is the basis for the SWA wind model discussed in 
Chapter III. Since the SAR integration time is constant, the distortion increases with the 
wind speed, eventually limiting the azimuthal resolution of the SAR at a finite cut-off 
wavenumber. This cut-off wavenumber can then be related to the wind speed. 

B. NORCSEX 95 

1. General 

NORCSEX ‘95 was conducted off the southwest coast of Norway over a twenty day 
period during September and October 1995. The third in the NORCSEX series, its’ 
objectives were to improve our understanding of the SAR imaging mechanisms of marine 
boundary layer processes such as fronts, eddies, and winds in coastal areas. Researchers from 
several organizations participated in the exercise including: Nansen Environmental and 
Remote Sensing Center, Bergen, Norway; Canada Center for Remote Sensing, Canada; Naval 
Postgraduate School, Monterey, USA; Environmental Research Institute of Michigan 
(ERIM), Ann Arbor, USA; and the French Research Institute for Exploitation of the Sea 
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(IFREMER), Brest, France. During the twenty day period there were a total of 14 satellite 
passes (7 per satellite) over the NORCSEX ‘95 area In-situ observations and data were 
collected on the R/V Hakon Mosby and a moored NPS buoy. The R/V Hakon Mosby 
attempted to stay in the SAR footprint and was used to collect data on features of particular 
interest in the SAR images in near real time. 

2. Platforms and Instruments 
a. R/V Hakon Mosby 

The University of Bergen operated the Research Vessel the R/V Hakon Mosby 
(Figure 2.7). It was instrumented and positioned to collect a variety of oceanographic and 
meteorological data within the SAR swath. The R/V Hakon Mosby left Bergen on 1 1 
September and arrived in the vicinity of 58 deg 30 min N, 4 deg 45 min E. The first 
operations were to deploy current meter and the NPS instrumented meteorology buoy. The 
orbits of both ERS-1 and ERS-2 during this period were such that the R/V Hakon Mosby 
could transit to all ascending and descending SAR swaths off the South and Southwest 
Norwegian coasts during the period. No swath was duplicated. 

The ocean surface and atmospheric conditions measured from the R/V Hakon 
Mosby that would be applicable to interpretation of SAR wind algorithms are SST, mean and 
turbulent wind, temperature and humidity, and mean pressure. The measured parameters on 
the R/V Hakon Mosby as well as on the NPS instumented buoy are listed in Table 2.1. All 
airflow measurements were made from a mast that was set up on a platform attached to the 
ship's main mast. The height of these sensors ranges from 15.3 to 18.7 m and were well 
exposed except for ship relative winds clockwise from 160 to 200 degrees. 
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Mean vector wind was measured with both a propeller-vane and a sonic 
anemometer, while mean temperature and humidity were obtained with naturally aspirated 
sensors. Ships course and speed (determined from both ship a gyro/log system and from GPS 
positions) were used to derive true vector wind from measured apparent winds. This 
information was needed because of frequent ship course and speed changes required from 
oceanographic surveying. Mean SST was measured by an infrared thermometer mounted 
above the starboard bridge. All mean quantities were sampled at 1 second rates and stored 
at 5 second averages. Turbulent wind and temperature were obtained from the sonic 
anemometer, sampled at 21 Hz. Turbulent humidity was measured with a IR 
shipboard fluctuating parameters were sampled at 21 Hz and stored 
processing for later analyses. 



hygrometer. All 
without further 
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Table 1: NORCSEX ‘95 Ship and Buoy 
Measurement Systems 




Ship 


Buoy 


Bulk Parameters (surface layer) 






Pressure 


X 


X 


Vector Wind 


X 


X 


Temperature (Air and Sea) 


X 


X 


Humidity 


X 


X 


Turbulence (surface layer) 






Vector Wind 


X 


X 


Temperature 


X 


X 


Humidity 


X 


X 


Surface Waves 






Directional Wave Spectra 




X 


Vertical structure (MABL) 






Vector Wind 


X 




Temperature 


X 




Relative Humidity 


X 





b. NPS Instrumented Buoy 

An instrumented moored meteorological and surface wave buoy, Figure 
2.6, was designed and constructed by NPS investigators especially for validation/ 
calibration of active (radar) remote sensors. The special features of the buoy were its 
capability to measure surface layer wind, wind stress, thermal stability and 2-dimensional 
surface wave properties from swell to centimeter wavelengths. The buoy hull was a toroid, 
2 m in diameter with a 5 m mast. The sensor and acquisition controls, mass memory 
storage, and batteries were housed in a cylindrical case (electronic case) mounted in the 
center of the toroid. The mooring chain and floatation design were sufficient to maintain 
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the buoy vertical under most wave conditions. 

Mean atmospheric surface layer properties were measured with a propeller- 
vane anemometer, and natural aspirated temperature and relative humidity sensors 
mounted near the 3 m level on the mast. Mean sea surface temperature was obtained from 
a thermocouple imbedded at the bottom of the buoy cylindrical electronics case. Mean 
pressure was measured by a sensor mounted in the electronics case. The mean data were 
sampled at a 1 Hz rate and averaged and stored at 5 -minute intervals. 

Turbulent wind and temperature were measured with a 3 -dimensional sonic 
anemometer mounted at the top of the mast, at the 5 m level. The surface wave 
conditions were measured with two separate systems. Wave conditions longer that the 2 
m toroid diameter, were interpreted from a system with three orthogonal accelerometers 
and three rate gyros. The resulting 3-dimensional motion information, sampled at 10 Hz, 
was also used to remove buoy motion contamination of sonic anemometer sensed 
velocities. Shorter waves, corresponding to capillary scale waves, were measured with 
three vertically oriented wires, spaced 4 cm apart, that penetrated the interface about 1 . 5 
meters upwind of the buoy. The turbulent and surface wave parameters were sampled at 
10 Hz rates and stored directly onto the mass memory devise for later time series analyses. 

The buoy was deployed from the R/V Hakon Mosby at 58 deg 30 min N 
and 4 deg 45 min E on 1 8 September. The mooring location was on the line between two 
SAR swaths, projected for the NORCSEX '95 period. It is believed that the buoy and its 
instrumentation systems worked properly for the 2-day period when the SAR swaths were 
near enough for its data to be used a surface truth. The complexity of the data analyses 



11 



limited its use for interpretations in this thesis. In particular, the merged buoy and sonic 
anemometer data analyses have not been developed enough to achieve full confidence in 
wind stress determination. 

c. ERS-2 

The ERS-2 satellite, launched in April 1995, is the second in a series 
designed primarily to provide high resolution information of the ocean environment. It is 
essentially the same as the ERS-1 with the exception of a number of enhancements and 
one new instrument designed to measure the chemical composition of the atmosphere 
(GOME). The other instruments carried by ERS-2 are the: 

- Active Microwave Instrument (AMI), 

- Radar Altimeter (RA), 

- Along Track Scanning Radiometer (ATSR), 

- Microwave Radiometer (MWR), and the 

- Precise Range and Range Rate Instrument (PRARE). 

The AMI incorporates the scatterometer and SAR and operates in three different modes to 
provide oceanic wind and wave measurements. The RA provides the precise altitude over 
all terrain. The ATSR measures sea surface temperatures and the ERS-2 version of this 
instrument has been enhanced by the addition of visible channels to monitor vegetation. 
The MWR and the PRARE are used to support the RA in determining the precise location 
of the satellite. 

Both ERS satellites are in a sun synchronous orbit crossing the equator on 
the descending node at about 10:30 local time each day. The satellites are currently in a 
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35 day repeat orbit with a one day separation with ERS-1 in the lead. This arrangement 
provides an excellent opportunity to compare the results from the two satellites. The 
orbits of either satellite can be changed if needed but it is expected that they will remain 
in this tandem mode. 

The instrument of interest in NORCSEX ‘95 is the AMI which 
incorporates the scatterometer and SAR and operates in three different modes to provide 
oceanic wind and wave measurements. It is the largest and uses the most power of all the 
instruments on board. The AMI has three modes of operation: scatterometer, SAR, and 
wave. The scatterometer mode can be used continuously and images a swath of 500 km 
on the sea surface. It provides wind speed and direction at a resolution of about 50 km 
and accuracies of 2 m/s and 20°. The SAR mode uses a C-band SAR to image a 100 km 
swath on the ocean surface. It consumes so much power that it can only be used for about 
ten minutes per orbit. The wave mode provides small images at 200 km intervals which 
are used to generate ocean-wave spectra, showing wave energy as a function of 
wavelength and direction. 
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Figure 2.1 Typical Satellite SAR Geometry (from Allan ed. 1983). 
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Figure 2.2 Bragg Scattering Geometry. Note the Most Preferential Location of 
Bragg Scatterers on the Larger Swell Wave Results in a Greater Backscatter 
Magnitude (as Indicated by the Magnitude of the Respective Reflected Signals) 
(after Curlander and Mcdonough, 1991) 
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Figure 2.3 Sketch Illustrating Both the Tilt and Hydrodynamic Mechanism. (From Alpers 
etal. 1981) 
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Figure 2.4 Modulation Transfer Function vs. Incidence Angle and Polarization, (from 
Alpers etal. 1981) 
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Figure 2.5 The Dependence of the Modulation Transfer Function 
Describing Tilt, Hydrodynamic, and Velocity Bunching Modulation on 
Azimuthal Angle. There is a Small Angular Interval Around the Range 
Direction Where Velocity Bunching is a Linear Process (From Alpers et 
al. 1981) 
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Figure 2.6 NPS Instrumented Buoy 





Figure 2.7 University of Bergen, Norway Research Vessel Hakon Mosby 
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III. WIND FIELD RETRIEVAL TECHNIQUES 



A. GENERAL 

Although there are many SAR wind algorithms, this thesis is focused on the two applied 
by Korsbakken (1996) on ERS-1 and ERS-2 SAR data from NORCSEX '95. Specifically, these 
were the CMOD4 wind algorithm, and the SAR Wind Algorithm (SWA). The description of the 
two algorithms following Korsbakken (1996) follows. 

In order to retrieve the absolute radar backscatter value (a„), necessaty to compute the 
wind field from the CMOD4 model, comprehensive calibration of the data is required, Scoon 
( 1 995). The calibration includes corrections for range spreading loss, for the antenna pattern, and 
for the power loss in analog to digital conversion (ADC). A brief description of these corrections 
is provided here to emphasize the importance of the calibration when using models which rely 

solely on the radar backscatter intensity. In extreme cases the values of o„ are corrected more 
than 5 dB. 

1. Range Spreading Loss 

The range spreading loss is the signal attenuation when the signal propagates from the 
radar to the target and back to the receiving antenna. The loss is expressed as: 

Q (dB ) = 1 01°g(-— — ) 

rL 



where R is the slant range distance to the target location and R ief is the reference slant range 
distance at mid swath (approximately 23°), Scoon (1995 ). The corrections for ERS-1 SAR 
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